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Abstract. We consider numerical approximation to the solution of non-autonomous 
evolution equations. The order of convergence of the simplest possible Magnus 
method will be investigated. 



1. Introduction 

The theoretical analysis of numerical methods and product formulae for partial 
differential equations plays an increasingly important role also in functional anal- 
ysis. In this paper we are interested in non-autonomous evolution equations of 
type 



(NCP S , X ) 



' ±u{t) = A{t)u{t), 
u(s) = x € X, 



t>s€ 



where X is a Banach space, (A(t), D(A(t))} is a family of (usually unbounded) linear 
operators on X. Such kind of problems arise in many applications in quantum 
physics, Hamiltonian dynamics, transport problems, etc. Among the numerical 
methods for the approximation of the solution, the Magnus integrators play an 
important role, see Isereles et al. [7] or Magnus [9]. The basic idea of this method 
is to express the solution u(t) in the form 

u(t) = exp(fl(t))x, 

where f2(t) is an infinite sum yielded by the formal Picard iteration 



[ A(r)dr~ [ \ f A(a)da,A(r 

JO ^ Jo [J 



(1) 



dr 



+ - 



A(n) d/j,,A(a) 



da, A(t) 



(It 



dr 



(with s = 0), where [U, V] = UV — VU denotes the commutator of the operators U 
and V. Numerical methods based on this expansion are presented in [8j by Iserles 
ct al. or in Hochbruck and Ostermann [5]. They are of the form 



(2) 



y n+1 = cxp(f2 n )y„, 
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to give an approximation of y(t n +x) at £ n +i = t n + h. Here fl n is a suitable 
approximation of tt(h) given by ([1]) when A{t) is substituted by + r). 

We want to investigate the convergence of the Magnus method where the mid- 
point rule is used, that is, in ^ 



In Section 2 we give a short overview on the theory of non-autonomous evolution 
equations. In Section 3 the norm convergence of the above Magnus method is 
proven, with an application to the Schrodingcr equation. Finally, in Section 4 we 
explore a slight gcncralizaton with applications to hyperbolic equations. 



In this section we summarize the main results and definitions on non- autonomous 
evolution equations and evolution semigroups needed for our later exposition. For 
a detailed account and bibliographic references see, e.g., the survey by Schnaubelt 
in [2j Section VI. 9.]. Consider now the non-autonomous Cuachy problem 



where X is a Banach space, (A(t), D(A(t))) is a family of (usually unbounded) 
linear operators on X . 

Definition 2.1. A continuous function u : [s,oo) — > X is called a (classical) 
solution of (NCP S>X ) if u e C x ([s, oo); X), u(t) £ D{A(t)) for all t > s, u(s) = x, 
and ^u(t) = A(t)u(t) for t > s. 

We use the following standard definition for the well-posedness of the non- 
autonomous Cauchy problem associated to a given family of linear operators. 

Definition 2.2. For a family (A(t), D(A(i))) tgffi of linear operators on the Ba- 
nach space X the non- autonomous Cauchy problem is well-posed (with regularity 
subspaces (y s ) se R and exponentially bounded solutions) if the following are true, 
(i) (Existence) For all s € M the subspacc 



Y s := < y £ X : there exists a classical solution for (NCP Sil ,) > C D(A(s)) 



is dense in X. 

(ii) (Uniqueness) For every y GY S the solution u s (-, y) is unique. 

(iii) (Continuous dependence) The solution depends continuously on s and y, i.e., 
if s n — > s € M, y n — > y £ Y s with y n £ Y Sn but in the topology of X, then we 
have 



(iv) (Exponential boundedness) There exist constants M > 1 and w£l such that 




2. Non-autonomous evolution equations 





\\u Sn (t,y n ) ~ u s (t,y)\\ 
uniformly for t in compact subsets of M, where 




|M*,y)||<Me^- s )||y|| 



for all y €Y S and t > s. 
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As in the autonomous case, the operator family solving a non-autonomous Cauchy 
problem enjoys certain algebraic properties. 

Definition 2.3. A family U = (U(t, s)) t > s of linear, bounded operators on a 
Banach space X is called an (exponentially bounded) evolution family if 

(i) U(t, r)U(r, s) = U(t, s), U(t, t) = I holds for all * > r > s G R, 

(ii) the mapping (t, s) n> U(t, s) is strongly continuous, 

(hi) ||?7(i 3 s)\\ < Me u ^- S ~> for some M > 1, w G R and alU > s G R. 

In general, however, and in contrast to the behavior of Co-semigroups (i.e., the 
autonomous case), the algebraic properties of an evolution family do not imply 
any differentiability on a dense subspace. So we have to add some differentiability 
assumptions in order to solve a non-autonomous Cauchy problem (abbreviated later 
on as NCP) by an evolution family. 

Definition 2.4. An evolution family U = (U(t, s))t> s is called evolution family 
solving NCP if for every s G R the regularity subspace 

Y s := jy G X : [s, oo) 3 t h-> U(t, s)y solves (NCP SjJ/ )} 

is dense in X. 

In this case, the unique classical solution of (NCP Sij; ) is given by u(t) := U(t, s)x. 
The well-poscdness of NCP can now be characterized by the existence of a solving 
evolution family, see [2j Proposition VI. 9. 3]. 

Proposition 2.5. Let X be a Banach space, and let (A(t),D(A(t))) be a family 
of linear operators on X. The following assertions are equivalent. 

(i) The associated non- autonomous Cauchy problem is well-posed, 
(ii) There exists a unique evolution family (U(t, s))t> s solving NCP. 

We say in this case that the evolution family is generated by the operators 
(A(t),D(A(t))) tm . 

Unfortunately, the well-poscdness of non-autonomous evolution equations is a 
complicated issue and there is no general theory describing it. Conditions implying 
well-poscdness arc generally divided into assumptions of "parabolic" and of "hyper- 
bolic" type. Roughly speaking, the main difference between these two types is that 
in the parabolic case we assume all A(t) being generators of analytic semigroups, 
while in the hyperbolic case we assume the stability for certain products instead. 
In both cases one has to add some continuity assumption on the mapping 1 1— > A(t). 
We mention only a typical and quite simple version for each type. 

Assumption 2.6 (Parabolic case). 

(PI) The domain D := D(A(t)) is dense in X and is independent of t G R. 

(P2) For each ( 6 1 the operator A{t) is the generator of an analytic semigroup 

e ,j4 W. For all t G R, the resolvent R(X,A(t)) exists for all A G C with 5RA > 

and there is a constant M > 1 such that 

||i2(A,A(t))||< |^ 

for 5RA > 0, t G R. The semigroups e ,A ® satisfy \\e sA ^\\ < Mc us for absolute 
constants u < and M > 1. 
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(P3) There exist constants L > and < a < 1 such that 

\\{A(t) - Ais^AiOy^l <L\t-s\ a for all t, s G M. 
Assumption 2.7 (Hyperbolic case). 

(HI) The family (A(i)) f6 R is stable, i.e., all operators are generators of Co- 
semigroups and there exist constants M > 1 and u£l such that 

(u>, oo) C p{A{t)) for all t G M 

and 

|]Ji?(A,A(^)) <M(A-w)- fc for all A > w 
i=i 

and every finite sequence — oo < t x < t 2 < • • • < < oo, fc G N. 
(H2) There exists a densely embedded subspace Y <^-t X, which is a core for every 
A(t) such that the family of the parts (A\ Y (t))t£«, in Y is a stable family on 
the space Y. 

(H3) The mapping K 9 t h-> A(t) £ C(Y, X) is uniformly continuous. 

Let us close this summary by recalling an important basic perturbation result, 
see Engel and Nagel [1 Theorem VI.9.19]. 

Theorem 2.8. Let {U(t,s))t> s be an evolution family in the Banach space X. Let 
B(t), t e R be closed operators such that U(t,s)X C D(B(t)), t i-> B(t)U(t,s) 
is strongly continuous and \\B(t)U(t, s) \\ < k(t — s) for t > s and some locally 
integrable function k. Then there is a unique evolution family (Usit, s))t> s such 
that 

(3) U B (t,s)x = U(t,s)x+ [ U B (t,r)B(r)U(r,s)xdr 



for all x G X and t > s. 

Note that this result only states the existence of an evolution family satisfying 
a variations-of constants formula, but not the well-posedness of non-autononous 
Cauchy problem associated to the operators (A(t) + B(t)) t ^K- 

3. Norm convergence of the Magnus method 

We would like to investigate the convergence of the following, simplest possible 
Magnus method introduced in the Introduction, where we take the first term of the 
Magnus series expansion and evaluate the integral using the midpoint rule. Let us 
take a time step h > 0, and consider the following iteration: 

(4) y„+i = e n "?/„, n n = hA (t n + f ) , 

starting with t = s, y = x and then t n+ \ = t n + h. Here and later on we will use 
the exponential notation to denote strongly continuous semigroups. 

Assumption 3.1. Assume that 

(a) there exists K > 1, u G M such that ||e tA ( r )|| < Ke ut for all r,i g 1. Further, 
A(t) = A + V(t), where A is the generator of a strongly continuous semigroup 
and V(t) G C(X) for all t G K, 

(b) ( Well-posedness) the non-autonomous Cauchy problem is well-posed, i.e., there 
is an evolution family U(t,s) which solves (NCP SjZ: ) and satisfies \\U(t, s)\\ < 
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(c) (Stability) the Magnus method j4| is stable, i.e., there are constants M > 1, 
uj G K such that for all n G N 



ft t- 



(2j + l)(t-s) ' 
In 



< Mc 



(fc+lWt-s) 



for all fc = 0, 1 . . . , n — 1, and 
(d) (Local Holder continuity) there exists a G (0, 1] and for all K > there is 
L = L(K) > such that 

\\A(t)-A(s)\\<L\t- S r 

for all t,se (-K,K). 
Here and later on, for bounded linear operators Lk G C(X), 



n-l 

n 

k=0 



L n —iL n -2 • ■ • Lq 



denotes the "time-ordered product". If the product is empty, we define it as the 
identity operator I on X. 

Note that the last condition (d) is always satisfied with a = 1 if we assume 
continuous differentiability of the map t M> V(t). Note also that the boundedness 
condition in (a) is always satisfied if the map 1 1— > V(t) is bounded. 

Theorem 3.2. Assume that the conditions in Assumption \3.1\ are satisfied. Then 
the Magnus method defined by (U) converges in the operator norm and has a uniform 
convergence order a G (0, 1] , i.e., for all t, s G M, t > s there is C > 1 such that 



U(t,s)- Je' 



lA(s+(2fe+l)t=a) 



< 



c 



for all n G N. 



Proof. First note that A(t) = A(s+±) + B(t) with B(t) = V(t)-V(s + ^) G C(X). 
Hence, we can apply the variation of constants formula ([3]) and obtain that 

U(s + h, s)x - e hA (s+h/2) x = I + u( ^ s + ]h r ) B (r) c (r-s)A(s+h/2) xAr 



Here the left-hand side is the difference of the exact solution and the solution 
obtained by the Magnus method after one time step (the so-called local error). 
Clearly, by the assumptions on V, we have 

||B(r)|| = \\V(r)-V(s + h/2)\\ < L\r - s - h/2\ a < Lh a 

for r G (s, s + h). Hence, 

(5) U(s + h,s)x~e hA{s+h ^x < J + Ke^+^L^Ke^-^WxWdr 

= K 2 Lc" h h a+1 \\x\\, 
proving the so-called consistency estimate we need for the following step. 
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The convergence order of the method can be established now using the standard 
telescopic argument. Let h := Then 

n-l 

U(t, «) = II U ( s + ( k + l ) h > s + kh ) 

k=0 

and the Magnus approximation is 

n-l 

W(t,s) = Y[c hA ^ +( - 2k+1 ^\ 

k=0 

Hence, for any x G X, 
U(t,s)x - W{t,s)x = 

n — 1 n—1 

= J2 U(s + (k + l)h,s + kh)\u(s + (j + l)h,s + jh)-e hA ^ +( - 2j+1) ^ 

j=0 k=j+l 

1=0 

From ([5]) we have 

U(s + (J + l)h, s + jh)x - e hA ( s+ ^ + ^^x < K 2 Le«> h h a+1 \\x\\, 
hence, for n sufficiently large, we obtain that 

\\U(t,s) - W(t,s)\\ <^2Ke {n - : >- 1)ujh ■ K 2 Le uh h a+1 ■ Me jujh 

3=0 

= K 3 LMh a+x ne nuh = K 3 LM(t - s )^ e ^-^) • — , 



since n = t-r^ . □ 

h 

Example 3.3. Motivated by Hochbruck and Lubich [3} Section 8.], we consider 
the d-dimensional Schrodinger equation 

= -\^ + b{x,t)il>, x£R d ,t>0, 
\tp(x,0) = ip (x), 

where the potential b and the initial function -0o are assumed to be 27r-periodic in 
every spatial coordinate. 
Assume further that 

• the function b is real-valued, 

• b(-,t) G i°°(M d ) for all t > 0, and 

• the functions &(■, t) G i°°(IR d ) is a- Holder continuous for some a G (0,1]. 
Then the operator family A(t) = |A — ib satisfies the conditions of Assumptions 
12.71 hence the Schrodinger equation is well-posed. Clearly, Assumptions 13 . 1 1 are also 
satisfied, hence the Magnus method Q has operator norm convergence of order a. 

Note that Hochbruck and Lubich showed that under additional smoothness con- 
ditions on 6, one can recover the classical second order convergence for smooth initial 
values. This convergence will be, however, in general only strong, the uniform rate 
remains first order. 
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We can use the above result to prove a regularity result on the propagators. Note 
that the compactness of propagators is important, for example, in case you want 
to establish a Peano type existence result on semilinear equations. 



Corollary 3.4. Assume that the conditions in Assumption \ 3.1\ are satisfied and 
that there is tq G R such that A(r$) generates an immediately compact semigroup. 
Then all the propagators U(t,s) are compact. 

Proof. We know by the bounded perturbation theorem that all the semigroups 
e tA ( r ) are immediately compact for r£l. Since the Magnus method ((4]) converges 
in the operator norm, it follows that the operators U (t, s) are compact. □ 

4. Generalization to strong convergence 

If we weaken the assumptions, the previous result can be generalized in a straight- 
forward way. We weaken the assumptions so that we only assume A(t) = A + V(t) 
where V(t) £ C(D(A), X) is allowed to be an unbounded operator. 

Assumption 4.1. Assume that 

(a) the non-autonomous Cauchy problem is well-posed and solved by the evolution 
family U(t, s), the operator A(r) is the generator of a strongly continuous semi- 
group for all r £ R, D := D(A(r)), all the graph norms arc equivalent with the 
same constants, and there is K > 1, u € R such that 

\\U(t,s)\\ <ife" (t - s) , ||c tA M|| <Ke ut 

for all t,rel, 

(b) (Stability) the Magnus method Q is stable, i.e., there are constants M > 1, 
to e R such that for all n G N 



k t -s ,( (2j + l)(t-s) 
A \ s+ 2n 



(k+l)uj(t-s) 

< Me n 



He v 

3=0 

for all k = 0, 1 . . . , n — 1, and 
(c) (Local Holder continuity) there exists a € (0, 1] and for all K > there is 
L = L(K) > such that 

\\A(t)-A(s)\\ D <L\t-s\ a 

for all t, s € (— K, K), where || ■ \\d denotes (one of) the graph norm(s) on D. 



Proposition 4.2. Assume that the conditions in Assumption \4-l\ are satisfied. 
Then the Magnus method defined by (j4]) converges for all x £ X and has for smooth 
initial values x G D convergence order a £ (0, 1]. 



The proof goes on the same line as that of Thcorem l3.2l and is therefore omitted. 
Example 4.3. Consider the Schrodingcr equation 



where, formally, 



if = -A(t)i>, x eR d ,t> 0, 
ip(x,0) = tpo(x), 



d 

A(t) : Yl d jaj , k (;t)d k 
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is the differential operator corresponding to the real- valued coefficients aj.k- See 
Engel and Nagel [2j Section VI. 5. c] for further details on such operators. We assume 
here that the coefficients aj.k satisfy 

a j , k (- > t) = a kd (- > t)€W 1 >°°(R d ), 

and 

d 

a j , k {x,t)y j y k > c\y\ 2 

j;k = l 

for a uniform constant c > and for all x,y € M. d . Then D(A(r)) = H 2 (M. d ) with 
graph norms equivalent to the Sobolcv norms. Differential operators of this type 
usually appear as the linearization of non-autonomous Schrodingcr equations. 

Assume further that there is a G (0,1] such that the maps t H> aj ik {-,t) are 
Q-Holder continuous. Since all the operators A(r) are selfadjoint, iA(r) generates 
a unitary group, hence the hyperbolicity Assumption 12 . 71 are satisfied. 

Hence, for all tp G H 2 (R d ) the Magnus method @J converges and its order is at 
least a. 
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